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A stability criterion is derived for self-similar solutions with perfect fluids which obey the equation 
of state P = kp in general relativity. A wide class of self-similar solutions turn out to be unstable 
against the so-called kink mode. The criterion is directly related to the classification of sonic 
points. The criterion gives a sufficient condition for instability of the solution. For a transonic point 
in collapse, all primary-direction nodal-point solutions are unstable, while all secondary-direction 
nodal-point solutions and saddle-point ones are stable against the kink mode. The situation is 
reversed in expansion. Applications are the following: the expanding flat Friedmann solution for 
1/3 < fc < 1 and the collapsing one for < fc < 1/3 are unstable; the static self-similar solution 
is unstable; nonanalytic self-similar collapse solutions are unstable; the Larson-Penston (attractor) 
solution is stable for this mode for < fc <■ 0.036, while it is unstable for 0.036 < fc; the Evans- 
Coleman (critical) solution is stable for this mode for < fc J; 0.89, while it is unstable for 0.89 ^ fc. 
The last application suggests that the Evans-Coleman solution for 0.89 fc is not critical because 
it has at least two unstable modes. 

PACS numbers: 0420D, 0425D, 0440, 0440N 



I. INTRODUCTION 



The invariance under scale transformation is one of the important features of gravitational physics as long-range 
force in general relativity as well as in Newtonian gravity. The scale invariance of field equations implies the existence of 
solutions which are invariant under the scale transformation. Such scale-free solutions are called self-similar solutions. 
Among them, the most widely researched system in general relativity is a spherically symmetric self-similar spacetime 
with a perfect fluid, which is pioneered by Cahill and Taub |^] . The existence of such solutions requires the equation 
of state to be the form P = kp \i it is barotropic. A wide class of matter fields can be modelled by this equation 
of state. In general relativity, spherically symmetric self-similar solutions are classified with the theory of dynamical 
systems [|| and with their asymptotic forms ||^. For a recent review, see [Q. 

Although self-similar solutions form a zero-measure set of solutions to the field equations, several authors conjectured 
that self-similar solutions play important roles in gravitational collapse and/or in cosmological situations (see, e.g., [^). 
In this context, the so-called Larson-Penston self-similar solution has been shown to describe generic gravitational 
collapse of isothermal gas in Newtonian gravity |^-^ . This solution was discovered by Penston and independently by 
Larson , and after that generalized to general relativity by Ori and Piran . Actually, Newtonian self-similar 

solutions are obtained by taking the limit fc ^ of those in general relativity pd| . Recently, it has been shown by 
numerical simulations and a normal mode analysis that the general relativistic counterpart of the Larson-Penston 
solution is an attractor solution in general relativistic gravitational collapse [l^ . This convergence phenomena, i.e., 
the generic convergence to a self-similar solution, can be regarded as the recovery of scale-invariance symmetry in 
gravitational physics because the scale-invariance symmetry which the system of the field equations originally has is 
regained in gravitational collapse even if the symmetry is lost in the initial situation. In the Larson-Penston solution, 
a naked singularity develops from regular initial data for < fc < 0.0105 ||l^,|l^. See also [|l2|,^ for the relationship 
between the cosmic censorship conjecture and the nature of the Larson-Penston solution as an attractor. 

In addition to the attractive nature of a self-similar solution, the critical behavior in the gravitational collapse 
discovered by Choptuik Q shed light on self-similarity (for a recent review, see ||l^). He investigated the threshold 
between the collapse to a black hole and the dispersion to infinity in the self-gravitating system of a massless scalar 
field. He found critical behavior, such as the mass-scaling law for the formed black hole mass, which is analogous 
to that in statistical physics. He also found that there is a discrete self-similar solution that sits at the threshold 
of the black hole formation, which is called a critical solution. After that, Evans and Coleman fl^ observed similar 
phenomena in the collapse of a radiation fluid (fc = 1/3) and found that the critical solution coincides with one of 
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the continuous self-similar solutions with analyticity, which we call Evans-Coleman solution [|. Recently, their work 
was extended to < A: < 1 IT^. A renormalization group approach by Koike et al. [|l8| gave a simple physical 
explanation to the critical phenomena and showed that the critical solution is characterized by a single unstable 
mode for a radiation fluid. After that their work was extended to < fc < 0.889 The detailed characteristics 

of the Evans-Coleman solution for < fc < 1 was investigated |l^j2^. Moreover, 'Hunter (a)' self-similar solution, 
which was originally discovered by Hunter was identified with the Newtonian counterpart of the Evans-Coleman 
solution §Jl|l. 

For Newtonian self-similar solutions, Ori and Piran [p2[ pointed out the existence of kink instability in a wide class 
of self-similar solutions. The kink mode results from the existence of sonic points. The nature of sonic points was 



examined in Newtonian gravity [p3| and in general relativity [|l l| , p4 - 26 1 . The kink mode corresponds to the injection 
of a weak discontinuity at sonic points at some initial moment. The instability is characterized by the divergence of 
the discontinuity in a finite time interval. The blow-up of this mode implies shock wave formation. Unfortunately, 
the analysis was restricted to Newtonian gravity and applies at most to the case < fc <C 1 In this paper, the 
analysis is generalized to fully general relativistic case. The stability criterion obtained here applies even to the case 
for any value of fc (0 < fc < 1). Then, we can deal with highly relativistic situations, such as the collapse of a radiation 
fluid, the evolution of early universe, and so on. The obtained criterion gives a sufficient condition for instability and 
a necessary condition for stability, because the analysis is specified on the kink mode. Of course, the present analysis 
reproduces the previously obtained results in Newtonian gravity. 

The organization of this paper is the following. In Section II, basic equations are presented. In Section III, the 
classification of sonic points are reviewed, which is closely related to the kink instability. In Section IV, the stability for 
the kink mode is analyzed in full order, and a stability criterion for this mode is derived. In Section V, applications to 
known self-similar solutions, such as the flat Friedmann solution, static self-similar solution, Larson-Penston solution 
and Evans-Coleman solution, are presented. In Section VI, we comment on recent numerical works by Neilsen and 
Choptuik and discuss two interesting limiting cases. Section VII is devoted to summary. We adopt units such 
that G = c=l. 

II. BASIC EQUATIONS 
A. Einstein's Equations 

In this section, we basically follow the formulation on a spherically symmetric self-similar system given by Bicknell 
and Hcnriksen |Q. It is noted that there are several different formulations (see, e.g., ) . 

The line element in a spherically symmetric spacetime is given by 

ds^ = -e^^'-^'Ut^ + e"(*''')dr2 + R^^t, r){d0'^ + sin^ 0d(j)^). (2.1) 

We consider a perfect fluid as a matter fleld 

T^"" = {p + P)u^'u'' + Pg"". (2.2) 

Here we adopt the comoving coordinates. Then Einstein's equations and the equations of motion for the matter field 
are reduced to the following partially differential equations (PDE's): 

da 2 dP 



dr p + P dr ' 

duj 2 dp A dR 



(2.3) 
(2.4) 



dt p + p at R at' 

'd^R ldaaR\ 1 ^da aP m , 

e"'^ — — = --^ -47rPi?, (2.7) 



dt^ 2 dt dt 2 ar ar R^ 



^ Although this self-similar solution itself had been already discovered by Ori and Piran we call this solution, together 
with its generalization to more general values of fc, Evans-Coleman solution for convenience in this paper. 
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and 



R 
"2 



dR 
'dt 



dR 
dr 



(2.8) 



where m{t,r) is called Misner-Sharp mass and five of the above six equations are independent. 
We assume the following equation of state: 

P = kp, 

where we assume < fc < 1. We define dimensionless functions such as 



V = 


Snr^ p, 


s = 


R 






r 




2m 


M = 






r 



We also define the following self-similar coordinate: 



It is often more convenient to use the following coordinates: 

r=-ln|i|, 
z = -ln|^|. 



(2.9) 

(2.10) 
(2.11) 

(2.12) 
(2.13) 



(2.14) 
(2.15) 



It is noted that ^ is a two-valued function of z, such as 2: = —00 corresponds to ^ = ±00 and z ~ +00 corresponds to 



(2.16) 
(2.17) 



It is found that equations ( p.3[) and ( |2.4[) are integrated as 



where aa{t) and aui{r) are arbitrary functions. These functions correspond to the freedom of re-scaling the time and 
radial coordinates as t = t{t) and f = r(t). Using this freedom, we let and a^^ be constant. 
For later convenience, we define a quantity y as 



y ^ 



M 



(2.18) 



which has a clear physical meaning that it is one third of the ratio of the 'averaged density' of the region interior 
to {t,r) to the local density at {t,r). We also define two velocity functions V and Vr. The V is the velocity of the 
^ = const curve relative to the fluid element, which is written as 



V = sign(^)e^+"2'' ^ 

while Vr is the velocity of the R = const curve relative to the fluid element, which is written as 

.S + S' 



Vr 



-e 2 



where 



(1) 

m 



d_ 

d_ 

dz 



V 



S + S' 



(2.19) 
(2.20) 

(2.21) 
(2.22) 
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V is written as 



y2 ^ a^^^^-^y2l^g- 

a„ 



Using these dimensionless variables, equations ( |2.5| )-(2.8) are reduced to the foUowing form: 

M + M' = -^kf]S^{S + S'), 



S + 2S' + S" 
k 







+ 1 




V J 





{S + S') 



1 + k 



52 



and 



M 
'J 



1 + e-^+^'-{S + S'f - e-^'iS + S'y, 



l\2 



(2.23) 



(2.24) 
(2.25) 



(2.26) 



(2.27) 



where three of the above four equations are independent. 

B. Self-similarity 

For self-similar solutions, we assume that all dimensionle ss qu a ntitie s depend only on z, i.e., M = M{z), S = S{z), 
rj = rj{z)^ a = <7{z), and uj = ll}{z). Then, from equations (2.24)-(2.27), we obtain the following ordinary differential 
equations (ODE's): 



S' 

v' 



1 + k y 

2(1 -y)- 2^^^ 



ky~ i(l + fc)2e'^r/ 



-k 



and constraint equation 

V^{1 - yf -{k + yf + (1 + kfe'^S-^{l - yf]S^) = 0. 
From equations ( ^.20 ) and ( ^.29[ ), we find 



Vb. = -V 



k + y' 



(2.28) 
(2.29) 
(2.30) 

(2.31) 

(2.32) 



III. SONIC POINTS 

The classification of sonic points in general relativity was studied by Bicknell and Henriksen ( [Q; see 
also [pl] , p5|j26| ) . Here we briefly review their work, which turn out to be closely related to the stability criterion 
of the kink mode. 



A. Regularity 



The sonic point is defined by 



V' = k, 



(3.1) 
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that is, V is equal to the sound speed. At the sonic point, the system of the ODE's (2.2S)-(2.30) is singular. Since 
V is the relative velocity of the ^ = const curve to the fluid element, informations with the sound speed propagate 
only outward or inward for V'^ > k, while in both directions for < fc, in terms of ^. In this sense, the former is 
supersonic, while the latter is subsonic, in terms of ^. Hence, informations can cross the sonic point only in a single 
direction. This is a physical reason why the system of the ODE's becomes singular at the sonic point. 
First, we consider the regularity condition at the sonic point. The regularity requires 



From equations ( ^.31 ), (3.1) and (3^), we find the following explicit relations at the sonic point 

{k + y,)^ + k{3ys - + Vs) 



' (1 + fc) i-fc ki (4ys) i-fc 



l + 3fc 



a-aCiuji'^ + k)^ 



g 2(l + fc) ' 



a„{l + kY 



l + fc 



(3.2) 



Therefore, sonic points are parametrized by one parameter y^- Equation (3.3) requires the following condition: 



Vs > yi 



'2k + ^fc(l + 6fc-3fc2) 
l + 3k ' 



(3.3) 

(3.4) 
(3.5) 

(3.6) 
(3.7) 



B. Classification of sonic points 



We introduce new independent variable u which is defined as 



du 



k. 



(3.8) 



Using u in place of 2, a sonic point turns out to be a stationary point (or a singular point) of the resultant system of 
the ODE's as 



dM 
du 

dS 1 — ?/ ,-,^2 



du 1 + k 
drj 
du 



l + fc y 

{V^ - k)S. 



iy' - k)M, 



2(1 - y){V^ -k)-2{ky- ^(1 + k^e^T^ 



V, 



A stationary point is classified by the behavior of solutions around the point. We put 

S = Ss{l + X2), 

e = 6(1 + 2:4), 



(3.9) 
(3.10) 
(3.11) 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



where 0:1 to X4 are regarded as components of the vector x. Then, we find the following linearized ODE's in the 
matrix form: 



—X = Ax, 
du 



(3.16) 
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where the matrix A is given by 



2k 

A = 



1 + k 



/ 2fc(l-2/) fc(l-fc)(l- v)/(l + fc) k{l-k){l~y)/{l + k) 

-2(1 -y) _(l-fc)(l-y)/(l + fc) -(l-fc)(l-y)/(l + fc) 

l + /c (4-3y)(l + /c) 2{1- k-y) 2(l-fc)(l-y) 

V -2(1 + -(l-'fc) -(1-^) 



(3.17) 



This matrix has two zero eigenvalues and two generically nonzero eigenvalues X± (A_ < X+) 



where X± are given by 



9k ~ 

A. = — i:A,. (3.18) 



with 

D = l + 6k-3k^ -iy- 24fcy + 12fc\ + iy^ + 12ky^. (3.20) 

These two eigenvalues A± are associated with the corresponding eigenvectors e±, respectively. The sonic point is 
classified into a saddle for A_ < < A+, a nondegenerate node for A_ < A_|_ < 0, a degenerate node for A_ = A+, 
and a focus for complex values of A±. For a saddle, only two solutions pass through it, one along the +ive direction 
e+ and the other along the — ive direction e_. For a nondegenerate node, there are one-parameter family of solutions 
which cross the sonic point along the -l-ive direction, while only one solution crosses along the —ive direction. The 
directions + and — for a node are called the primary and secondary directions, respectively. For a degenerate node, 
two directions + and — coincide with each other. The sonic point which is a focus is unphysical. The sonic point is 
a nondegenerate node or saddle for £> > 0, a degenerate node for D = 0, and a focal point for D < 0. The condition 
D > requires 

Vs < y2 or 2/3 < Vs, (3.21) 

where 2/2,3 (2/2 < ys) are given by 

1 + 6A; - 3fc2 ip ^3fc(l-fc)(l + 6fc-3fc2) 



and D = holds only for ys = j/2 or y^ = y^. Together with equation (3.7), since yi < y2 can be verified, it is found 
that 

yi <ys < y2 or ya < ys, (3.23) 

for the sonic point to be physical. Whether the sonic point is a node or saddle is determined by the sign of [(1+k)'^ — D]. 
It changes its sign at y^ = y^ or y^ = y5 (y4 > ys), where 



1 + 6A: - 3fc2 ± V(l + 4fc-F)2 + 8fc2(l-fc)2 

y^^^ = 2(i+3fc) — ^ — - 

Since ys < yi can be verified, it is found that the sonic point is unphysical for ys < yi or y2 < ys < ys, a degenerate 
node for ys = y2 or ys = ya, a nondegenerate node for yi < ys < y2 or ys < jjs < yi, and a saddle for < ys- The 
classification of sonic points in terms of the parameter ys is plotted in figure H^. 
It is seen that, along the allowed directions ±, 

^ . ^~k-2y + X^ ^3^^^ 
77 1 — fe 

at the sonic point. Using this and the relation 

V 1-kii 3 - fc - 2y 



V 1 + k r] 



(3.26) 
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which follows from equations (2.23) and (2.29), it is found that 



{V^y = -A±, (3.27) 



V _i Vd 

y-2^2(TTfc)' ^^-^^^ 

where and hereafter the upper and lower signs denote the +ive and — ive directions, respectively. Therefore, if 
a solution crosses a saddle sonic point from subsonic interior to supersonic exterior, then it goes along the —ive 
direction. A sonic point with subsonic interior and supersonic exterior is called a 'transonic' point. If we consider 
a solution with regular center and sonic points, the solution must have at least one transonic point. In principle, a 
solution may have other kinds of sonic points, such as the one with supersonic interior and subsonic exterior which we 
call an anti-transonic point. It is found that the sonic point with supersonic interior and subsonic exterior must be a 
saddle and the solution passes it along the -l-ive direction. Furthermore, it is found that V' /V < 1/2 or {V"^)' < k for 
the -|-ive direction, V' /V > 1/2 or {V"^)' > k for the —ive direction, and V /V — 1/2 or {V"^)' = A: for a degenerate 
node. 



C. Analyticity 



Furthermore we can require the analyticity condition (Taylor-series expandability) at the sonic point, i.e. 

M = Ms -f M,,i(z - z,) + Ms,2{z - zsf • • • , 

S ^ Ss+ SsAz - Zs) + Ss,2{z ~ Zsf H , 

rj = i]s + r]s.i{z - Zs) + ijsM^ ~ ^sf H , 



around the sonic point. From equations (2.28)-(2.31), we find 



Ms 1 = — ^M., 

'YTk^'' 

_ {3k + 4ys-5)TVD 
"^^'^ ~ 2(1 - fc) ''^ 



(3.29) 
(3.30) 
(3.31) 



(3.32) 
(3.33) 
(3.34) 



if D is not negative. It is found that the upper and lower signs correspond to the -|-ive and —ive directions, respectively. 
In other words, there are two analytic solutions at a sonic point, one of which crosses the sonic point along the -|-ive 
direction and the other along the —ive direction. Since all coefficients of Taylor series are determined only by M^, Ss 
and r\s^ the analytic solution with each direction is unique for given j/j. Hence, for a saddle, two solutions pass it and 
both of them are analytic. On the other hand, for a nondegenerate node, one analytic solution and one-parameter 
family of nonanalytic solutions pass it along the primary direction, while one solution passes it along the secondary 
direction and this solution is analytic. In particular, all nonanalytic solutions must cross nodes and go along the 
primary direction. 



IV. KINK INSTABILITY 
A. Equations for kink mode 

We consider perturbations which satisfy the following conditions in the background of a self-similar solution: (1) 
The initial perturbations vanish inside the sonic point for t < 0. (For t > 0, the initial perturbations vanish outside 
the sonic point.) (2) M, S and are continuous everywhere, in particular at the sonic point. (3) 77' is discontinuous 
at the sonic point, although it has definite one-sided values as z — ^ — and z ^ Zs + 0. 

We denote the full order perturbations as 
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SS{t, z) = S{t, z) - Sb{z), 
dM{T,z) = M{t, z) - Mb{z), 
St]{t,z) = ?7(t,z) - ?7b(z), 



(4.1) 
(4.2) 
(4.3) 



where Mb, Sb and rjb denote the background self-similar solution. Hereafter we often omit the s ubscript '6' unless the 
omission may cause confusion. By conditions (2) and (3), and from equations ( ^.24 ) and ( ^.25 ), it is found that SM' 
and 6S' are also continuous. Hence, the perturbations satisfy 5M = 0,SS = 0, 6r] = 0, 6M' = 0, 6S' — 0, Srj' 7^ at 
the sonic point at initial moment r = tq, where the prime denotes the derivative with respect to z on the perturbed 
side. The evolution of the initially unperturbed region is completely described by the background self-similar solution 
because no information from the perturbed side can penetrate the unperturbed side by condition (1). Then, by 
conditions (2) and (3), we find SM = 0, (55 = 0,6r] = 0,6M = 0,SS = O^Sfj = 0,SM' = 0,(55" = 0,(5?)' 0, at the 
sonic point for r > tq for the case oi t < (for r < tq for the case oi t > 0). This mode has a physical meaning 
that it injects discontinuity into the density gradient. In fact, the density gradient with respect to the physical length 
dl = e'^^'^dr on the t = const spacelike hypersurface is directly related to 77' as 



dp 
dl 



dr 



V 



1 



(4.4) 



Differentiating equations ( p. 24 ) and ( 2.25D with respect to z and estimating both sides at the point of discontinuity, 
we obtain the following equations for the full order perturbations: 



6M" 
M 
SS" 
S 



k I5i 



l + ky rj 

y ^v' 

1 + k rj ' 



(4.5) 
(4.6) 



We also differentiate equation (2.26) with respect to z. The result is 




5^ 
5 



+ k7]' + 3fc?7— 



We find from equations ( 2.26 ) and (4/7) for the full order perturbations at the point of discontinuity 

k 5' Srj' kV-^ 



6S" 

SS" 



1 + k S r/ 
(55"' k 



1 + k 

kv-^ 

1 + k 



2k 
1 + k 



5 1 + k 
V 

l-kSi]' 
1 + k rj 

t 
V 



1 + k 

Sr]' St]" 

V V 
SS" 

5 ^ 



(4.7) 



(4.8) 



T]' 6r,' 



(5?7' 5" k 



'S 

5i 6S" 



1 + k rj S 1 + k rj 5 



1 



1 + k jj 



S 

'I 



1 



1 + fc 
St]' 



1 



5 



if' 



1 



5 
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(5r/" ^77' Si]' Srj'^ 



S 



h — 
77 



S 



1] 



- 2 



S' 

6S" 
S 



6S" 



At the sonic point, equation (18) is reduced to 

S 



y Si]' 

1 + k rj 



= 0, 



(4.9) 



(4.10) 



which is equivalent to equation (4.(:). Therefore, it is found that the injection of this type of discontinuity into the 
density gradient field i s po ssible only at the sonic point. Next, we evaluate both sides of equation (^^) at the sonic 
point using equations (3.1) and (3.2). One can observe that there appear SS'" and Sr]" in this equation. It is easily 
found that such terms appear in the equation only through the following combination: 



SS" 



s ■ 1 + . , • ^'-''^ 

In order to write these terms using the lower order derivatives, we eliminate M' from equations ( p. 24 ) and (2.25) as 



M^M ^ rjS^ + (1 + k)rjS^S' + krjS^S. 
Differentiating the above equation twice with respect to z, we have 
M" - M" 

+ (1 + k){f]"S^S' + 4r]'SS'^ + 2rj'S^S" + 2r)S'^ + QrjSS'S" + r^S^S'") 
+k{'n"S^S + 4f]' SS'S + 2r]'S^S' + 2r]S'^S + 2'qSS"S + Aj]SS'S 

Estimating the full order perturbations of both sides at the sonic point, we find 



(4.12) 



TjS^S") 



(4.13) 



SS 


/// 


y 


Stj" 


s " 


i + k 


V 


y 


'SM" 


SM' 


1 ^ 


^k 


M 


M 




3 




(- + 




1 + k 





k SS" 
l + k~ 



SS" 

s 



- 2 



3 S' 
1 + fc^ 

St]' SS" 
rj S 



S^ 
S 



Sjj' 
V 



(4.14) 



Using the above relation, we evaluate both sides of equation (4.9) at the sonic point using SM", SS" and Srj'. After 
a rather lengthy calculation, we obtain the following equation: 



SJf 



-2 



l-krj' b - ik - Ay\ Si l-k f Srj' 
'1 + k rj ^ 1 + k J rj 1 + 



— I =0. 



(4.15) 



It is noted that the PDE's become an ODE for the discontinuity in the density gradient at the sonic point. This is 
due t o the nature of the sonic point. It is more convenient to rewrite this in terms of V and SV' using the relation 
( 3.26 ). The discontinuity in V' is directly related to the discontinuity in Srj' through 

SV _ l-k Srj' 

'V ^ "i + fc V 

Finally we obtain the following equation for the full order perturbation as 



(4.16) 



SV 

'V 



1 



SV 

~v~ 



SV 

~v~ 



0. 



(4.17) 



It is noted that the above equation has complete correspondence to its Newtonian counterpart, equation (15) of p2[ , 
when the latter is rewritten in terms of the relative velocity of the constant self-similar coordinate curve to the fluid 
element. 
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B. Linear order analysis 



By linearizing equation (4.17), we can obtain 



5V' 



l-2(^ 
V 



This equation is integrated as 



where 



const • e 



'V 

a-l-2(- 



(4.18) 



(4.19) 



(4.20) 



Therefore, it is found that the discontinuity in V' grows for a > 0, decays for a < and the mode is neutral for a — 
as T increases. On the other hand, it is found that the discontinuity in V' grows for a < 0, decays for a > and the 
mode is neutral for a = as t decreases. The above discussion will suggest some kind of instability. Nevertheless, 
the linear order analysis will not be sufficient in this case from the following reason. The discontinuity in the density 
gradient which is injected as a perturbation may not be included in the background solution. Then, there is no natural 
measure how much the perturbation grows compared with the background solution. In the linear order analysis, we 
can show at most a power-law growth of the discontinuity for sufficiently small perturbation. In the next subsection, 
we will show a divergence of full order perturbation at some finite moment, even if initial perturbation is small. This 
does indicate instability. 



C. Full order analysis 



In fact, equation ( 4.17 ) is easily integrated. Putting 



X 



5V' 



we rewrite equation ( 4.17 ) as 



X - + = 0. 

There are two stationary solutions X ~Q and X — a. General solutions are the following: 

a 



X 



X = 



1 — exp(— ar + const 

■, for a = 



, for a 7^ 0, 



T + const 

For an initial value X = Xq 7^ at t = tq, the solution is written as 



1 _ 1 



-a(T— ro) 



1 1 

a Xr 



, for a 7^ 0, 



^ = ^ + (r - To), for a = 0. 



(4.21) 
(4.22) 

(4.23) 
(4.24) 

(4.25) 
(4.26) 
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1. D > case 



First, we consider a > 0. From equation (4.25), it is seen that, as r increases, X blows up to — oo at some finite 
moment r ^ — > tq for Xq < 0, wliere 



1 , / 1 a 

= To + - In 1 - — 

a \ Xq 



(4.27) 



while X monotonically approaches a for Xq > 0. Next, we consider a < 0. From equation ( 4.25 ), it is seen that, as r 
increases, X monotonically approaches for Xq > a, while X blows up to —oo at some finite moment t = — > tq 
for Xq < a. 

Actually these two analyses under two different backgrounds are only apparently different local pictures of nonlinear 
dynamics. If we consider the background solution which crosses the sonic point along the +ive direction, a is positive 
and the stationary solution X = a corresponds to another possible direction, since 



V 



The situation is reversed for the — ive direction solution, since 

1 - 2 



(4.28) 



(4.29) 



The global picture is the following. With the initial value V /V = {V' /V)q at t = tq, we obtain the following full 
order behavior of the perturbation: 



V 

YL 

V 
V 



—oo as T —r Td — for 
const for 

oo for 



as T 



< 

I 

V 



< 



as T increases, while 



v 




V 




v 




V 




v 






— > oo 


V 





as T 



-oo for 



< 



'V 

const for ( — 





V 

< 



(4.30) 
(4.31) 
(4.32) 

(4.33) 
(4.34) 
(4.35) 



as T decreases. The dynamical behavior of the perturbation is schematically depicted in figure 



2. D = case 



We consider a = 0. From equation (4.26), it is seen that, as r increases, X diverges to — oo at some finite moment 
T ^ Td, where 



Td = TQ 



1 

'x'q'- 



(4.36) 



for Xq < 0, while X approaches for Xq > 0. 

This case corresponds to a degenerate node for which D = Q and {V' /V)b = 1/2. Then we find 
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as T increases, while 



V fV'\ 1 

— ^-oo as T^Td-O for { — \ < -, (4.37) 

V I 1 fV'\ 
- as T^oo for - < ( — (4.38) 



V 2 2 - \ V 



V 1 fV'\ 1 

> - as r ^ -c» for — < -, (4.39) 

v 1 /y\ 

— ^oo as + for 2<(l7) ' ^"^'^^^ 



as T decreases. The behavior is depicted in figure ||. 



D. Stability criterion 



First, we consider t < Q. The time evolution is regarded as the increase of t from r = tq. Then, we obtain the 
following criterion: all +ive direction solutions and degenerate-nodal-point ones are unstable, while all — direction 
solutions are stable for this mode. Next, we consider the case t > 0. In that case, the time evolution is regarded as the 
decrease of t from t = tq. Then, we find the following criterion: all — direction solutions and degenerate-nodal-point 
ones are unstable, while all -I— direction solutions are stable for this mode. These are summarized in table |[ 

The stability criterion which applies to all possible situations is the following: if V' >V/2 at the sonic point, then 
the solution is unstable, while, if not, then the solution is stable for the kink mode. This criterion can be rewritten 
in terms of the partial derivative with respect to i or r as follows: in terms of the derivative with respect to if 
dV^ /dt-\- k/t > 0, then the solution is unstable, while, if not, the solution is stable for the kink mode; in terms of the 
partial derivative with respect to r, if dV/dr > V/{2r) at the sonic point, the solution is unstable, while, if not, the 
solution is stable for this mode. 

However, it should be noted that we cannot say that a solution is stable because our analysis is specified on the 
kink mode. The present criterion for stability should be considered as a necessary condition for stability, while the 
criterion for instability should be considered as a sufficient condition for instability. 



V. APPLICATIONS 



A. Flat Friedmann solution 



The fiat Friedmann solution is a member of self-similar solutions. This solution is characterized by y = 1/ 3 and 



ry'/r/ = 2. This solution crosses a transonic point. For j/^ = 1/3, we find two nonzero eigenvalues from equation ( 3.1£ ) 
as 

A = 4-^4. (5.1) 

Therefore, the sonic point is a nondegenerate node for /c 7^ 1/3 and a degenerate node for k — 1/3. Since we find from 
equation { \i-2l ) 



A = -A:-i, (5.2) 

the fiat Friedmann solution crosses a nondegenerate node along the primary direction forO< k < 1/3, a degenerate 
node for fc = 1/3, and a nondegenerate node along the secondary direction for 1/3 < A: < 1. The solution is expanding 
for t > 0, while the solution is collapsing for t < 0. From the obtained stability criterion, we find that the expanding 
fiat Friedmann solution is unstable against the kink mode for 1/3 < fc < 1, which implies that it may not be a 
good cosmological model for 1/3 < k < 1. We find that the collapsing fiat Friedmann solution is also unstable for 
< A; < 1/3 against the kink mode, which implies that generic gravitational collapse does not proceed homogeneously 
for < fc < 1/3. We also find that the expanding fiat Friedmann solution for < A: < 1/3 and the collapsing fiat 
Friedmann solution for 1/3 < fc < 1 do not suffer kink instability. 
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B. Static self-similar solution 



The static self-similar solution is characterized hy y — I and t]' /t] — 0. This solution crosses a transonic point. For 



Us = 1, we find eigenvalues from equation (3.19): 

A = -2fc,-l + fc. (5.3) 

Since we find from equation ( ^.25 ) 

A = -l + fc, (5.4) 

the static self-similar solution crosses a nondegenerate node along the secondary direction for < fc < 1/3, a degenerate 
node for k — 1/3, and a nondegenerate node along the primary direction for 1/3 < /c < 1. From the stability criterion 
for the kink mode, we find that the solution is unstable for i > for < /c < 1/3. We also find that the solution is 
unstable for i < for 1/3 < fc < 1. For fc = 1/3, the solution is unstable, irrespective of the sign of t. However, it 
should be noted that the moment t — has no physical meaning in the static solution because the system is invariant 
under time translation. Therefore, we conclude that the static self-similar solution is unstable for < fc < 1. 

C. Nonanalytic self-similar solutions 

As we have seen, there are one-parameter family of nonanalytic self-similar solutions which cross a sonic point. All 
of them cross nodes along the primary direction. When the stability criterion is applied, these nonanalytic solutions 
are all unstable for t < 0. Thereby, stable self-similar collapse solutions with sonic points must be analytic. If 
we consider self-similar solutions with regular center and sonic points, a set of all stable collapse solutions becomes 
discrete, while a set of all collapse solutions is dense. In other words, the stability requirement for the kink mode 
considerably reduces the number of self-similar collapse solutions. For t > 0, nonanalytic self-similar solutions do not 
suffer kink instability. 

D. Anti-transonic point 

We have seen that an anti-transonic point must be a saddle and that the solution crosses it along the -|-ive direction. 
From the stability criterion, the solution is unstable for t < 0. Therefore, in collapse, any stable self-similar solution 
cannot pass an anti-transonic point. 

E. Larson-Penston (attractor) solution 

The Larson-Penston solution is a solution which is analytic both at the center and at a transonic point and it has no 
zero in the velocity field Vr. The Larson-Penston solution has no analytic unstable mode in Newtonian gravity 
and in general relativity at least for < fc < 0.03 [|l^. The nature of the sonic point was studied by |]ll] , po| . 

This solution crosses the sonic point along the secondary direction for < fc ^ 0.036. It has been found that 
the sequence of these analytic solutions changes from the secondary-direction nodal-point solution to the degenerate- 
nodal-point one at fc ~ 0.036. For a larger value of fc, the sequence changes from the degenerate-nodal-point solution 
to the primary-direction nodal-point one|^ By the present stability criterion for i < 0, we conclude that the Larson- 
Penston solution is stable for the kink mode for < fc ^ 0.036, while it is unstable for 0.036 ^ fc. It is important that 
the kink mode does not affect the nature of the Larson-Penston solution as an attractor for < fc < 0.036, while it 
does for 0.036 < fc. 

For t > 0, the time-reversed Larson-Penston solution is unstable against the kink mode for < fc < 0.036, while it is 
stable for 0.036 < fc for this mode. This suggests that the solution cannot describe a realistic expanding inhomogeneous 
universe for < fc < 0.036. From the Larson-Penston solution for t < 0, there is another method to construct the 
solution for t > 0. That is the analytic continuation of the solution for t < to t > 0. Since the Larson-Penston 
solution is 'quasi-static', the analytic continuation beyond ^ = is possible [p HlO|JTl| , p6t . This continuated solution for 
t > does not cross a sonic point |p^ , ^ , ^ . Therefore, the solution has no kink mode for t > 0. 



^ Although we could not find that in fiterature, we have confirmed it numerically. 
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F. Evans-Coleman (critical) solution 



Evans-Coleman solution is a solution which is analytic both at the center and a transonic point and it has one 
zero in the velocity field Vr. The Evans-Coleman solution has a single analytic unstable mode, which was shown for 
< fc < 0.889 by normal mode analyses |]l8|,^ and for < fc < 1 by numerical simulations The nature of a 

sonic point and the asymptotic form of this solution were studied by 0,|2^ ■ 

The sequence of these analytic solutions changes its character for k ~ 0.41 and k ~ 0.89 in terms of the sonic point. 
The solution crosses a saddle sonic point for < fc ^ 0.41. For 0.41 ^ fc ^ 0.89, the solution crosses a nodal sonic 
point in the secondary direction, a degenerate-nodal sonic point for fc ~ 0.89, and a nodal sonic point in the primary 
direction for 0.89 ^ fc. From the stability criterion, it is found that the kink mode does not affect the critical nature 
of the Evans-Coleman solution for < fc ^ 0.89, while it is also found that the Evans-Coleman solution suffers the 
kink instability for 0.89 ^ fc. Since the critical solution is assumed to have a single unstable mode, it indicates that 
the Evans-Coleman solution for 0.89 ^ fc cannot be a critical solution because the solution has one analytic unstable 
mode and one non-analytic unstable mode, i.e., the kink mode. 

For t > 0, the time-reversed Evans-Coleman solution is unstable for < fc ^ 0.89, while it is stable for 0.89 ^ fc for 
the kink mode. This suggests that the solution would not be a good model of the expanding inhomogeneous universe 
for < fc < 0.89. The behavior of the solution as ^ increases from — oo changes at fc ~ 0.28 from the quasi-static 



solution to the 'asymptotically Minkowski' solution as fc increases |20 . Only for the quasi-static solution, the analytic 
continuation beyond f = is possible. However, it was shown that the analytically continuated Evans-Colcman 
solution beyond ^ = encounters another sonic point and will not be able to cross it regularly [ ^J20|j2^ . 



VI. DISCUSSIONS 



A. Comment on Neilsen and Choptuik |17| 



Here we should comment on the recent work by Neilsen and Choptuik ^Tj. They numerically simulated the PDE's 
(Einstein's equations and equations of motion for a perfect fluid) and found critical behavior even for 0.89 fc- 
They found that the critical solution obtained by solving the PDE's well agreed with the Evans-Coleman self-similar 
solution which is obtained by solving the ODE's with the requirement of analyticity at the sonic point. By numerical 
simulations of the PDE's, they measured the critical exponent of the power law which the formed black hole mass 
obeys. The obtained value of critical exponent is continuous with respect to the change of fc. It implies the existence 
of one analytic unstable mode. In contrast, we have shown that the Evans-Coleman solution suffers kink instability for 
0.89 < fc. We would like to point out the following three possibilities. The first is that perturbations which belonged 
to the kink mode in their numerical simulations might have been so small that there would have been little time for 
the kink mode to grow up sufficiently until the numerical simulations were stopped for other reasons. The second 
is that some kind of numerical viscosity, which is usually included in numerical codes of fluid dynamics deliberately 
or not deliberately, might have killed the kink instability. The third possibility is that the nonlinear coupling of the 
kink mode with other analytic modes might have weakened the growth of instability. In spite of such possibilities, for 
0.89 ^ fc, the Evans-Coleman has two unstable modes, one is analytic and the other is not analytic. 

B. Two Limiting cases 

It is remarkable that the stability analysis for the kink mode of self-similar solutions in general relativity is completely 
parallel to that in Newtonian gravity. The stability for the kink mode is determined by the class of the pertinent 
sonic point. Therefore, it is clear that the present analysis on the stability of general relativistic self-similar solutions 
contains the Newtonian analysis for isothermal gas as a limit of fc — > 0. On the other hand, it should be noted that 
fc = 0, i.e., a dust fluid, is not included in the present analysis since self-similar solutions with a dust fluid cannot be 
regarded as a continuous limit of fc ^ of those with P = fcp in several respects ||lT| . 

Since we have concentrated on self-similar solutions with P — kn for < fc < 1, w e have not dealt with a 
'stiff matter', i.e., fc = 1. This is because, as seen in equations (3^)-(3^), ( 3.25 ), ( 3.26| ), ( |3.34 ) and ( 4.16| ), it is 



clear that more careful treatment is needed for that case. The stiff matter is one of the most important examples 
that can be described by the perfect fluid model because it is equivalent to a massless scalar field under certain 
conditions pTf . In smte of this equivalence, the critical phenomena which have been observed in these two systems 
look very different [[l4|jl7|]. Hence, the analysis for a stiff matter will be very interesting. 
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VII. SUMMARY 



We have studied the stabihty of self-similar solutions with perfect fluids in general relativity. It has been found 
that a wide class of self-similar solutions turn out to be unstable against the kink mode. The development of this 
instability will result in the formation of a shock wave. Since the present analysis is specified only on the kink mode, 
we cannot say about generic outcome of the growth of this instability. However, it is probable that the nonlinear 
coupling with analytic modes may play important roles in the development of the instability. The proposition that 
the stability of solutions for this mode is solely governed by the nature of the pertinent sonic point applies to general 
relativity as well as to Newtonian gravity. 

The obtained criterion has been applied to several known self-similar solutions in general relativity. Then, it 
turns out that the expanding flat Friedmann solution has a growing kink mode for 1/3 < fc < 1. The existence of 
the growing kink mode in the flat Friedmann solution with a radiation fluid might affect the standard cosmological 
structure formation scenario. We have also investigated the stability of self-similar collapse solutions which have 
recently called attention in the convergence and critical phenomena in the gravitational collapse of a perfect fluid in 
general relativity. The results are that the Larson-Penston (attractor) solution loses its attractive nature for 0.036 < k 
and that the Evans-Coleman (critical) solution loses its critical nature for 0.89 < k. 
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TABLE I. Stability for the kink mode and the class of sonic points 



t 


Primary node 


Secondary node 


Degenerate node 


Saddle: + + 


Saddle: - 


> 


Stable 


Unstable 


Unstable 


Stable 


Unstable 


< 


Unstable 


Stable 


Unstable 


Unstable 


Stable 



t No transonic point belongs to this class, while all anti-transonic points belong to this class. 
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FIG. 1. Classification of sonic points in terms of a parameter ye is depicted. The sonic point is unphysical for i/s < yi, a 
nondegenerate node for yi < ys < 2/2 or j/3 < j/s < 2/4, a focus (unphysical) for 2/2 < J/s < 2/3, a degenerate node for 2/s = 2/2 
or 2/s = y-i, and a saddle for 2/4 < ys- The flat Priedmann solution (2/s = 1/3) and static self-similar solution (2/s = 1) are also 
plotted. See text for details. 
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FIG. 3. Dynamical behavior of the kink mode for a degenerate-nodal sonic point. 
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